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Abstract 



The Y production cross section in proton-proton collisions at a/s = 7 TeV is measured 
using a data sample collected with the CMS detector at the LHC, corresponding to 
an integrated luminosity of 3.1 ± 0.3 pb^^. Integrated over the rapidity range \y\ <2, 
we find the product of the Y(1S) production cross section and branching fraction to 
dimuons to be cr{pp Y(IS)X) ■ 6(Y(1S) ^ ^^^') = 7.37 ± 0.13+°;^^ ± 0.81 nb, 
where the first uncertainty is statistical, the second is systematic, and the third is as- 
sociated with the estimation of the integrated luminosity of the data sample. This 
cross section is obtained assuming unpolarized Y(1S) production. If the Y(1S) pro- 
duction polarization is fully transverse or fully longitudinal the cross section changes 
by about 20%. We also report the measurement of the Y(1S), Y(2S), and Y(3S) dif- 
ferential cross sections as a function of transverse momentum and rapidity. 



Submitted to Physical Review D 



*See Appendix A for the list of collaboration members 



1 



1 Introduction 

The hadroproduction of quarkonia is not understood. None of the existing theories success- 
fully reproduces both the differential cross section and the polarization measurements of the 
]/tp or Y states [Ij. It is expected that studying quarkonium hadroproduction at higher center- 
of-mass energies and over a wider rapidity range will facilitate significant improvements in our 
understanding. Measurements of the Y resonances are particularly important since the theo- 
retical calculations are more robust than for the charmonium family, due to the heavy bottom 
quark and the absence of b-hadron feed-down. Measurements of quarkonium hadroproduction 
cross sections and production polarizations made at the Large Hadron Collider (LHC) will al- 
low important tests of several alternative theoretical approaches. These include non-relativistic 
QCD (NRQCD) factorization [2J, where quarkonium production includes color-octet compo- 
nents, and calculations made in the color-singlet model including next-to-leading order (NLO) 
corrections [3| which reproduce the differential cross sections measured at the Tevatron exper- 
iments mini without requiring a significant color-octet contribution. 

This paper presents the first measurement of the Y(1S), Y(2S), and Y(3S) production cross 
sections in proton-proton collisions at ^ys = 7 TeV, using data recorded by the Compact Muon 
Solenoid (CMS) experiment between April and September 2010. In these measurements the 
signal efficiencies are determined with data. Consequently, Monte Carlo simulation is used 
only in the evaluation of the geometric and kinematic acceptances. The document is organized 
as follows. Section |2] contains a short description of the CMS detector. Section |3] presents the 
data collection, the online and offline event selections, the Y reconstruction, and the Monte 
Carlo simulation. The detector acceptance and efficiencies to reconstruct Y events in CMS are 
discussed in Sections |4]and|5| respectively. In Section|6]the fitting technique employed to extract 
the cross section is presented. The evaluation of systematic uncertainties on the measurements 
is described in Section|7] Section |8] presents the Y(nS) cross section results and comparisons to 
existing measurements at lower collision energies ||4l|5| and to the pythia [6 J event generator. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m inner diameter, 
providing a field of 3.8 T. Inside the solenoid in order of increasing distance from the interac- 
tion point are the silicon pixel and strip tracker, the crystal electromagnetic calorimeter, and 
the brass/ scintillator hadron calorimeter. Muons are detected by three types of gas-ionization 
detectors embedded in the steel return yoke: drift tubes (DT), cathode strip chambers (CSC), 
and resistive plate chambers (RFC). The muon measurement covers the pseudorapidity range 
\fj\ < 2.4, where ?/ = — ln[tan(0/2)] and the polar angle 6 is measured from the z-axis, which 
points along the counterclockwise beam direction. The silicon tracker is composed of pixel 
detectors (three barrel layers and two forward disks on either side of the detector, made of 
66 million 100 ■ 150 }im^ pixels), followed by microstrip detectors (ten barrel layers plus three 
inner disks and nine forward disks on either side of the detector, with 10 million strips of pitch 
between 80 and 184 fim). Due to the strong magnetic field and the fine granularity of the sili- 
con tracker, the transverse momentum, pj, of the muons matched to silicon tracks is measured 
with a resolution of about 1% for a typical muon in this analysis. The silicon tracker also pro- 
vides the primary vertex position with ~ 20 f/m accuracy. The first level (LI) of the CMS trigger 
system, composed of custom hardware processors, uses information from the calorimeters and 
muon detectors to select the most interesting events. The high-level trigger (HLT) further de- 
creases the event rate before data storage. A more detailed description of the CMS detector can 
be found elsewhere HZj- 
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3 Data sample and event reconstruction 



3 Data sample and event reconstruction 
3.1 Event selection 

The data sample used in this analysis was recorded by the CMS detector in ppcoUisions at a 
center-of-mass energy of 7TeV. The sample corresponds to a total integrated luminosity of 

3.1 ± 0.3 pb^^ ||8|. The maximum instantaneous luminosity was 1031 cm s and event pileup 
was negligible. Data are included in the analysis if the silicon tracker, the muon detectors, and 
the trigger were performing well and the luminosity measurement was available. 

The trigger requires the detection of two muons at the hardware level, without any further 
selection at the HLT The coincidence of two muon signals without an explicit pj requirement 
is sufficient to maintain the dimuon trigger without prescaling. All three muon systems - DT, 
CSC, and RFC - take part in the trigger decision. 

Anomalous events arising from beam-gas interactions or beam scraping in the beam transport 
system near the interaction point, which produce a large number of hits in the pixel detector, are 
removed with offline software filters [9J. A good primary vertex is also required, as defined in 
Ref. lU. The detector systems are aligned and calibrated using LHC collision data and cosmic- 
ray muons l|TO| . 

3.2 Monte Carlo simulation 

Simulated events are used to tune the selection criteria and to compare with data. Upsilon 
events are produced using pythia 6.412 lU, which generates events based on the leading- 
order color-singlet and octet mechanisms, with NRQCD matrix elements tuned by comparing 
calculations with the CDF data [llj and applying the normalization and wavefunctions as rec- 
ommended in Ref. |T2|. The generation of Y(2S) and Y(3S) states has been included by chang- 
ing the relevant particle masses and branching ratios. The simulation includes the generation 
of Xb states. Final-state radiation (FSR) is implemented using PHOTOS lfT3llT4| . The response 
of the CMS detector is simulated with a GEANT 4-based flSl Monte Carlo (MC) program. The 
simulated events are processed with the same reconstruction algorithms as used for data. 

3.3 Offline muon reconstruction 

In this analysis, a muon is defined as a track reconstructed in the silicon tracker and associated 
with a compatible signal in the muon detectors. Tracks are reconstructed using a Kalman filter 
technique which starts from hits in the pixel system and extrapolates outward to the silicon 
strip tracker. The tracks found in the silicon tracker are propagated to the muon system and 
required to be matched to at least one muon segment in one muon station. Further details may 
be found in Ref. il6| . 

Quality criteria are applied to the tracks to reject muons from kaon and pion decays. The tracks 
are required to have at least twelve hits in the silicon tracker, at least one of which must be in 
the pixel detector, and a track-fit X'^ per degree of freedom smaller than five. In addition the 
tracks are required to emanate from a cylinder of radius 2 mm and length 50 cm centered on 
the pp interaction region and parallel to the beam line. The muons are required to satisfy: 

pl^ > 3.5GeV/c if |f/''| < 1.6 , or pt^ > 2.5GeV/c if 1.6 < |?/''| < 2.4 . (1) 

These kinematic criteria are chosen to ensure that the trigger and muon reconstruction efficien- 
cies are high and not rapidly changing within the acceptance window of the analysis. 



3.4 Y event selection 
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The momentum measurement of charged tracks in the CMS detector is affected by system- 
atic uncertainties caused by imperfect knowledge of the magnetic field, the amount of ma- 
terial, and sub-detector misalignments, as well as by biases in the algorithms which fit the 
track trajectory. A mismeasurement of the track momentum results in a shift and broadening 
of the reconstructed peaks of dimuon resonances. An improved understanding of the CMS 
magnetic field, detector alignment, and material budget was obtained from cosmic-ray muon 
and LHC collision data llTOl [TtI ITsU , and the residual effects are determined by studying the 
dependence of the reconstructed ]/xp dimuon invariant-mass distribution on the muon kine- 
matics ||T9| . The transverse momentum corrected for the residual scale distortion is parame- 
terized as Pt = {1 + ui + ci2^^) • Pj, where pj is the measured muon transverse momentum, 
fli = (3.8 ± 1.9) ■ 10^"*, and fl2 = (3.0 ± 0.7) ■ 10^"^. The coefficients for terms linear in rj and 
quadratic in pj and pj ■ rj are consistent with zero and are not included. 

3.4 Y event selection 

To identify events containing an Y decay, muons with opposite charges are paired, and the 
invariant mass of the muon pair is required to be between 8 and 14GeV/c^. The longitudinal 
separation between the two muons at their points of closest approach to the beam axis is re- 
quired to be less than 2 cm. The two muon helices are fit with a common vertex constraint, 
and events are retained if the fit probability is larger than 0.1%. The dimuon candidate is 
confirmed offline to have passed the trigger requirements. If multiple dimuon candidates are 
found in the same event, the candidate with the best vertex quality is retained; the fraction of 
signal candidates rejected by this requirement is about 0.2%. Finally, the rapidity, y, of the Y 
candidates is required to satisfy < 2 because the acceptance diminishes rapidly at larger ra- 
pidity. The rapidity is defined as i/ = ^ In ^ ^~^\ ) ' ^^^^^ ^ energy and p || the momentum 
parallel to the beam axis of the muon pair. 

The dimuon invariant-mass spectrum in the Y(nS) region for the dimuon transverse momen- 
tum interval pj < 30GeV/c is shown in Fig.[l]for the pseudorapidity intervals |f/''| < 2.4 (left) 
and \ f]f'\ < 1.0 (right). We obtain a Y(1S) mass resolution of 96 ± 2MeV/c^ including muons 
from the entire pseudorapidity range, and 69 ± 2MeV/c^ when both muons satisfy < 1. 
The observed resolutions, determined as a parameter of the fit function described in Section|6| 
are in good agreement with the predictions from MC simulation. 



4 Acceptance 

The Y — > acceptance of the CMS detector is defined as the product of two terms. The 

first is the fraction of upsilons of given pj and y such that each of the two muons satisfies 
Eq. ([T|. The second is the probability that each muon can be reconstructed in the tracker using 
the CMS software, in the absence of further event activity in the detector and without quality 
cuts imposed. Both components are evaluated by simulation and parametrized as a function 
of the Pj and rapidity of the Y. 

The acceptance is calculated from the ratio 

where Njg^ {pT,y) is the number of upsilons generated in a {pT,y) bin, while N^^, {pT,y) is the 
number reconstructed in the same {pT,y) region but now using the reconstructed, rather than 
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4 Acceptance 
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Figure 1: The dimuon invariant-mass distribution in the vicinity of the Y(nS) resonances for 
the full rapidity covered by the analysis (left) and for the subset of events where the pseudo- 
rapidity of each muon satisfies |f/''| < 1 (right). The solid line shows the result of a fit to the 
invariant-mass distribution before accounting for acceptance and efficiency, with the dashed 
line denoting the background component. Details of the fit are described in Section |6] 

generated, variables. In addition, the numerator requires that the two muons reconstructed in 
the silicon tracker satisfy Eq. ([T]). 

The acceptance is evaluated with a signal MC sample in which the Y decay to two muons is 
generated with the evtgen [20J package including the effects of final-state radiation. There 
are no particles in the event besides the Y, its daughter muons, and final-state radiation. The 
upsilons are generated uniformly in pj and rapidity. This sample is then fully simulated and 
reconstructed with the CMS detector simulation software to assess the effects of multiple scat- 
tering and finite resolution of the detector. The acceptance is calculated for two-dimensional 
(2-D) bins of size (1 GeV/c ■ 0.1) in the reconstructed pj and y of the Y and used in candidate- 
by-candidate yield corrections. 

The 2-D acceptance map for unpolarized Y(1S) is shown in the left plot of Fig. |2] The accep- 
tance varies with the resonance mass. This is shown in the right plot of Fig. |2| which displays 
the acceptance integrated over the rapidity range as a function of pj for each upsilon resonance. 
The transverse-momentum threshold for muon detection, especially in the forward region, is 
small compared to the upsilon mass. Therefore, when the Y decays at rest, both muons are 
likely to reach the muon detector. When the Y has a modest boost, the probability is greater 
that one muon will be below the muon detection threshold and the acceptance decreases until 
the Y transverse momentum reaches about 5 GeV/ c, after which the acceptance rises slowly. 
The production polarization of the Y strongly influences the muon angular distributions and 
is expected to change as a function of pj. In order to account for this, the acceptance is calcu- 
lated for five extreme polarization scenarios pT| : unpolarized and polarized longitudinally and 
transversely with respect to two different reference frames. The first is the helicity frame (HX), 
defined by the flight direction of the Y in the center-of-mass system of the colliding beams. The 
second is the CoUins-Soper (CS) frame Ii22j , defined by the bisection of the incoming proton 
directions in the Y rest frame. 
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Figure 2: (Left) Unpolarized Y(1S) acceptance as a function of pj and y; (Right) the unpolarized 
Y(1S), Y(2S), and Y(3S) acceptances integrated over rapidity as a function of pj. 

5 Efficiency 

We factor the total muon efficiency into three conditional terms, 

e(total) = e(trig|id) ■ e(id|track) ■ e (track [accepted) = etrig • £id • £track ■ (3) 

The tracking efficiency, £track/ combines the efficiency that the accepted track of a muon from 
the Y(nS) decay is reconstructed in the presence of other activity in the silicon tracker, as deter- 
mined with a track-embedding technique Il23l , and the efficiency for the track to satisfy qual- 
ity criteria, as determined with the tag-and-probe (T&P) technique l|23l described below. The 
muon identification efficiency, £id, is the probability that the track in the silicon tracker is iden- 
tified as a muon. It has been computed as described in Ref. [16.1 and is also based on the T&P 
method. The efficiency that an identified muon satisfies the trigger, etrig, is again measured 
with the same technique. 

The tag and probe technique is a data-based method used in this analysis to determine the 
track quality, muon trigger, and muon identification efficiencies. It utilizes dimuons from ]/xp 
decays to provide a sample of probe objects. A well-identified muon, the tag, is combined with 
a second object in the event, the probe, and the invariant mass is computed. The tag-probe pairs 
are divided into two samples, depending on whether the probe satisfies or not the criteria for 
the efficiency being evaluated. The two tag-probe mass distributions contain a ]/xp peak. The 
integral of the peak is the number of probes that satisfy or fail to satisfy the imposed criteria. 
The efficiency parameter is extracted from a simultaneous unbinned maximum-likelihood fit 
to both mass distributions. 

The ]/xp resonance is utilized for T&P efficiency measurements as it provides a large-yield and 
statistically-independent dimuon sample l|24| . To avoid trigger bias, events containing a tag 
and probe pair have been collected with triggers that do not impose requirements on the probe 
from the detector subsystem related to the efficiency measurement. For the track-quality effi- 
ciency measurement, the trigger requires two muons at LI in the muon system without using 
the silicon tracker. For the muon-identification and trigger efficiencies, the trigger requires a 
muon at the HLT, that is matched to the tag, paired with a silicon track of opposite sign and the 
invariant mass of the pair is required to be in the vicinity of the ]/xp mass. 

The component of the tracking efficiency measured with the track-embedding technique is well 
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Figure 3: Single-muon identification efficiencies as a function of pj for six \r]^'\ regions, mea- 
sured from data using J/t/^ T&P (closed circles). The efficiencies determined with Y MC truth 
(triangles), ]/xp MC truth (open circles), and ]/ip MC T&P (squares), used in the evaluation of 
systematic uncertainties, are also shown. 
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Table 1: Single-muon identification efficiencies, in percent, measured from J/j/^ data with T&P. 
The statistical uncertainties in the least significant digits are given in parentheses; uncertain- 
ties less than 0.05 are denoted by 0. For asymmetric uncertainties the positive uncertainty is 
reported first. 



Pt I'?^ 



(GeV/c) 


0.0-0.4 


0.4-0.8 


0.8-1.2 


1.2-1.6 


1.6-2.0 


2.0-2.4 


2.5-3.0 










100(0,4) 


94(6) 


3.0-3.5 










95(3) 


100(0,4 


3.5^.0 


83(2) 


89(2) 


88(2) 


96(3) 


100(0,2) 


100(0,4 


4.0-i.5 


92(2) 


95(2) 


99(1,3) 


98(2,3) 


100(0,2) 


100(0,6 


4.5-5.0 


99(1,2) 


99(1,3) 


95(3) 


96(3) 


100(0,2) 


100(0,4 


5.0-6.0 


98(2) 


100(0,1) 


100(0,2) 


100(0,1) 


97(3) 


100(0,5 


6.0-8.0 


100(0,2) 


100(0,1) 


100(0,1) 


100(0,2) 


100(0,2) 


94(6,7 


8.0-50.0 


100(0,2) 


97(3) 


100(0,3) 


97(3,4) 


100(0,3) 


98(2,9 



Table 2: Single-muon trigger efficiencies, in percent, measured from ]/ip data with T&P. The 
statistical uncertainties in the least significant digits are given in parentheses; uncertainties less 
than 0.05 are denoted by 0. For asymmetric uncertainties the positive uncertainty is reported 
first. 



Pt I'/'' 



(GeV/c) 


0.0-0.4 


0.4-0.8 


0.8-1.2 


1.2-1.6 


1.6-2.0 


2.0-2.4 


2.5-3.0 










93(1) 


92(2) 


3.0-3.5 










94(1) 


93(1) 


3.5-4.0 


69(1) 


81(1) 
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98(1) 


94(1) 


97(1) 


4.0^.5 


79(1) 
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92(1) 


96(1) 


4.5-5.0 


85(1) 


95(1) 
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97(1) 


96(1) 


99(1) 


5.0-6.0 


90(1) 
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85(1) 


99(0,1) 
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96(1) 


6.0-8.0 
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97(1) 


85(1) 
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8.0-50.0 


92(1) 


97(1) 


86(1) 


99(1) 


97(1) 


99(2) 



described by a constant value of (99.64 ± 0.05)%. The efficiency of the track-quality criteria 
measured by the T&P method is likewise nearly uniform and has an average value of (98.66 ± 
0.05)%. Tracks satisfying the quality criteria are the probes for the muon identification study. 
The resulting single-muon identification efficiencies as a function of pj for six \?]^'\ regions are 
shown in Fig.|3]and Table [l] The probes that satisfy the muon identification criteria are in turn 
the probes for the study of the trigger efficiency. The resulting trigger efficiencies for the same 
Pj and |?/^'| regions are shown in Fig. |4] and Table |2j 

Figures |3] and |4] also show single-muon identification and trigger efficiencies, respectively, de- 
termined from a high-statistics MC simulation. The single-muon efficiencies determined with 
the T&P technique in the data are found to be consistent, within the uncertainties and over 
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Figure 4: Single-muon trigger efficiericies as a function of Pj for six \ f]t^ \ regions, measured from 
data using ^/tp T&P (closed circles). The efficiencies determined with Y MC truth (triangles), 
]/tp MC truth (open circles), and ]/tp MC T&P (squares), used in the evaluation of systematic 
uncertainties, are also shown. 



most of the kinematic range of interest, with the efficiencies obtained from the Y MC simu- 
lation utilizing the generator-level particle information ("MC truth"). Two exceptions are the 
single-muon trigger efficiency for the intervals |//''| < 0.4 and 0.8 < \rj^\ < 1.2, where the 
efficiency is lower in data than in the MC simulation. In both cases the MC simulation is 
known to not fully reproduce the detector properties or performance: gaps in the DT coverage 
< 0.4) and suboptimal timing synchronization between the overlapping CSC and DT 
subsystems (0.8 < < 1.2). For all cases the data-determined efficiencies are used to obtain 
the central results. 

The Y efficiency is estimated from the product of single-muon efficiencies. Differences between 
the single and dimuon efficiencies determined from MC truth and those measured with the 
T&P technique can arise from the kinematic distributions of the probes and from bin averaging. 
This is evaluated by comparing the single-muon and dimuon efficiencies as determined using 
the T&P method in ]/xp — ?► MC events to the efficiencies obtained in the same events 

utilizing generator-level particle information. In addition, effects arising from differences in 
the kinematic distributions between the Y and ]/^p decay muons are investigated by comparing 
the efficiencies determined from Y — >• ji^fi^ MC events to those from ]/xp — > fi^fi^ MC events. 
In all cases the differences in the efficiency values are not significant, and are used only as an 
estimate of the associated systematic uncertainties. 

The efficiency of the vertex probability cut is determined using the high-statistics ]/xp data 
sample, to which the Y selection criteria are applied. The efficiency is extracted from a simul- 
taneous fit to the dimuon mass distribution of the passing and failing candidates. It is found to 
be (99.2 ± 0.1)%. A possible difference between the efficiency of the vertex probability cut 
for the ]/ip and Y is evaluated by applying the same technique to large MC signal samples of 
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each resonance. No significant difference in the efficiencies is found. The efficiency of the re- 
maining selection criteria listed in Section|3]is studied in data and MC simulation and is found 
to be consistent with unity. 



6 Measurement of the cross sections 

The Y(nS) differential cross section is determined from the signal yield, N^*, obtained directly 
from a weighted fit to the dimuon invariant-mass spectrum, after correcting for the acceptance 
(A) and the total efficiency (e), through the equation 

'P-iPP^^^-Sm.s^YinS)^,^,-) (4) 



dpj dy 



C -ApT-Ai/' 



(5) 



upon normalization by the integrated luminosity of the dataset, C, and by the bin widths, Apj 
and Ay, of the Y transverse momentum and rapidity. 

The Y(1S), Y(2S), and Y(3S) yields are extracted via an extended unbinned maximum likeli- 
hood fit. The measured mass-lineshape of each Y state is parameterized by a "Crystal Ball" 
(CB) function [25 1; this is a Gaussian resolution function with the low side tail replaced with a 
power law describing FSR. The resolution, given by the Gaussian standard deviation, is a free 
parameter in the fit but is constrained to scale with the ratios of the resonance masses. The 
FSR tail is fixed to the MC shape. Since the three resonances overlap in the measured dimuon 
mass, we fit the three Y(nS) states simultaneously. Therefore, the probability distribution func- 
tion (PDF) describing the signal consists of three CB functions. The mass of the Y(1S) is a free 
parameter in the fit, to accommodate a possible bias in the momentum scale calibration. The 
number of free parameters is reduced by fixing the Y(2S) and Y(3S) mass differences, rela- 
tive to the Y(1S), to their world average values |I26|. A second-order polynomial is chosen to 
describe the background in the 8-14GeV/c^ mass-fit range. 

The fit to the dimuon invariant-mass spectrum, before accounting for acceptance and efficien- 
cies, is shown in Fig. [l]for the Y transverse momentum interval pj < 30GeV/c, and for the 
fifteen pj intervals used for the Y(1S) differential cross-section measurement in Fig. [s] The 
observed Y(nS) signal yields are reported in Table |3] The width of the pj intervals chosen for 
each resonance reflects the corresponding available signal statistics. In all cases the quoted un- 
certainty is statistical. As shown in Table |3j for each resonance the sum of the yields in each 
Pt interval is consistent with the yield determined from a fit to the entire px range. Given the 
significant rj and px dependencies of the efficiencies and acceptances of the muons from Y(nS) 
decays, we correct for them on a candidate-by-candidate basis before performing the mass fit 
to obtain N^'' in Eq. (jiji. Specifically: an Y candidate reconstructed with px and y from muons 
with pj^'^ and //''i-^ is corrected with a weight 

ZV = Wacc ■ 2^track ' ^id " ^trig " ^misc (6) 

where the factors are: (i) acceptance, zvacc = 1/ A^ {pj,y)', (ii) tracking, Wfrack = 1/ ^track' 

(iii) 

identification, t^id = 1/ [£id{pj\r]^^) ■ eid(Px'/ '/^'')] ; M trigger, Wtng = 1/ [etrig(Px ' etrig(Px' 
and (v) additional selection criteria, Wmisc/ including the efficiency of the vertex selection crite- 
ria. The acceptance depends on the resonance mass; the Y(3S) gives rise to higher-momenta 
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6 Measurement of the cross sections 
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Figure 5: Fit to the dimuon invariant-mass distribution in the specified px regions for \y\ < 
2, before accounting for acceptance and efficiency. The solid line shows the result of the fit 
described in the text, with the dashed line representing the background component. 
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muons which results in a roughly 10% larger acceptance for the Y(3S) than for the Y(1S). Con- 
sequently, the corrected yield for each of the Y(nS) resonances is obtained from a fit in which 
the corresponding Y(nS) acceptance is employed. 

We determine the Y(nS) differential cross section using the above procedure separately for 
each polarization scenario. The results are summarized in Table |4] We also divide the data 
into two ranges of rapidity, \y\ < 1 and 1 < \y\ < 2, and repeat the fits to obtain the Y(nS) 
differential cross sections reported in Table |5] The integrated cross section for each resonance 
is obtained from the corresponding sum of the differential cross sections. The results for the 
Y(1S) px-integrated, rapidity-differential cross section are shown in Table|6] 

7 Systematic uncertainties 

Systematic uncertainties are described in this section, together with the methods used in their 
determination. We give a representative value for each uncertainty in parentheses. 

We determine the cross section using acceptance maps corresponding to five different polar- 
ization scenarios, expected to represent extreme cases. The values of the cross section obtained 
vary by about 20%. The variations depend on pj thus affecting the shapes of the pj spectrum. 

The statistical uncertainties on the acceptance and efficiencies - single-muon trigger and muon 
ID, quality criteria, tracking and vertex quality - give rise to systematic uncertainties in the 
cross-section measurement. We vary the dimuon event weights in the fit coherently by ± lt7"(stat. ) 
The muon identification and trigger efficiencies are varied coherently when estimating the as- 
sociated systematic uncertainties (8%). 

The selection criteria requiring the muons to be consistent with emanating from the same pri- 
mary vertex are fully efficient. This has been confirmed in data and simulation. The selection 
of one candidate per event using the largest vertex probability also has an efficiency consistent 
with unity. We assign an uncertainty (0.2%) from the frequency of occurrence of signal can- 
didates in the data that are rejected by the largest vertex probability requirement but pass all 
the remaining selection criteria. The muon charge misassignment is estimated to be less than 
0.01% |27) and contributes a negligible uncertainty. 

Final-state radiation is incorporated into the simulation using the PHOTOS algorithm. To es- 
timate the systematic uncertainty associated with this procedure, the acceptance is calculated 
without FSR and 20% of the difference is taken as the uncertainty based on a study in Ref . ||T4| 
(0.8%). 

The definition of acceptance used in this analysis requires that the muons from the Y decay 
produce reconstructible tracks. The kinematic selection is applied to the reconstructed pj and 
rj values of these tracks. Uncertainties on the measurement of track parameters also affect the 
acceptance as a systematic uncertainty. The dominant uncertainty is associated with the mea- 
surement of the track transverse momentum. The acceptance is sensitive to biases in track mo- 
mentum and to differences in resolution between the simulated and measured distributions. 
The magnitude of these effects is quantified by comparing measurements of resonance mass 
and width between simulation and data I.19J . To determine the effect on the Y acceptance, we 
introduce a track pj bias of 0.2%, chosen to be four times the maximum momentum scale resid- 
ual bias after calibration (0.3%). We also vary the transverse momentum resolution by ±10%, 
corresponding to the uncertainty in the resolution measurement using ]/xp, and recalculate the 
acceptance map (0.1%). 
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8 Results and discussion 



Imperfect knowledge of the production pj spectrum of the Y resonances at ^ys = 7 TeV con- 
tributes a systematic uncertainty. The Y MC sample used for the acceptance calculation, Eq. 
was generated flat in pj, whereas the pj spectrum in the data peaks at a few GeV/ c, and be- 
haves as a power law above 5 GeV/c. To study the effect of this difference, we have re-weighted 
the sample in pj to more closely describe the expected distribution in data based on a fit to the 
spectrum obtained from PYTHIA (1%). 

The distribution of the z position of the ppinteraction point influences the acceptance. We have 
produced MC samples of Y(nS) at different positions along the beam line, between —10 and 
+10 cm with respect to the center of the nominal collision region (1%). 

High-statistics MC simulations were performed to compare T&P single-muon and dimuon ef- 
ficiencies to the actual MC values for both the Y and ]/tp, see Figs. |3] and |4] The differences and 
their associated uncertainties are taken as a source of systematic uncertainty. The contributions 
are: possible bias in the T&P technique (0.1%), differences in the ]/xp and Y kinematics (1%), and 
taking the product of single-muon efficiencies as an estimate of the double-muon Y efficiencies 
(1.6%). 

MC trials of the fitter demonstrate that it is consistent with providing an unbiased estimate 
of the yield of each resonance, its mass, and the mass resolution (1%). A systematic variation 
may arise from differences between the dimuon invariant-mass distribution in the data and 
in the PDFs chosen for the signal and background components in the fit. We consider the 
following variations in the signal PDF. As the CB parameters which describe the radiative tail 
of each resonance are fixed from MC simulation in the nominal fit to the data, we vary the 
CB parameters by three times their uncertainties (3%). We also remove the resonance mass 
difference constraint in the pj integrated fit (0.6%). We vary the background PDF by replacing 
the polynomial by a linear function, while restricting the fit to the mass range 8-12 GeV/c^ (3% 
when fitting the full px and y ranges, varying with differential interval). 

The determination of the integrated luminosity normalization is made with an uncertainty of 
11% [8J. The relative systematic uncertainties from each source are summarized in Table [7| 
for the full rapidity range, for two rapidity ranges in Table |8| and for five rapidity ranges in 
Table |9] The largest sources of systematic uncertainty arise from the statistical precision of the 
T&P determination of the efficiencies from the data and from the luminosity normalization 
which dominates. 

8 Results and discussion 

The analysis of the collision data acquired by the CMS experiment at = 7 TeV, correspond- 
ing to an integrated luminosity of 3.1 ± 0.3 pb^^, yields a measurement of the Y(nS) integrated 
production cross sections for the range \y\ < 2: 

(7(pp ^ Y(IS)X) • S(Y(1S) ^ fi+fi-) = 7.37 ±0. 13 (stat.)+°;^^(syst.) ±0.81 (lumi.)nb, 
(7(pp ^ Y(2S)X) ■ ;B(Y(2S) ^ = 1.90±0.09(stat.);°;^0(syst.) ±0.24(lumi.)nb, (7) 

(7(pp ^ Y(3S)X) ■ B(Y(3S) ^ = 1.02±0.07(stat.);g;J^(syst.) ±0.11(lumi.)nb. (8) 

The Y(1S) and Y(2S) measurements include feed-down from higher-mass states, such as the Xb 
family and the Y(3S). These measurements assume unpolarized Y(nS) production. Assump- 
tions of fully-transverse or fully-longitudinal polarizations change the cross sections by about 
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20%. The px-differential Y(nS) cross sections for the rapidity intervals |t/| < 1, 1 < < 2, and 
lyl < 2 are shown in Fig.|6j The px dependence of the cross section in the two exclusive rapid- 
ity intervals is the same within the uncertainties. The Y(1S) pj-integrated, rapidity-differential 
cross sections are shown in the left plot of Fig. [7| The cross section shows a slight decline to- 
wards ly I =2, consistent with the expectation from pythia. The ratios of Y(nS) cross sections 
differential in px are reported in Table 10 and shown in the right plot of Fig.[7[ The uncertainty 
associated with the luminosity determination cancels in the computation of the ratios. Both 
ratios increase with px. In Fig.|8]the differential cross sections for the Y(1S), Y(2S), and Y(3S) 
are compared to PYTHIA. The normalized px-spectrum prediction from PYTHIA is consistent 
with the measurements, while the integrated cross section is overestimated by about a factor 
of two. We have not included parameter uncertainties in the PYTHIA calculation. We do not 
compare our measurements to other models as no published predictions exist at = 7 TeV 
for Y production. 
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Figure 6: Y(nS) differential cross sections in the rapidity interval \y\ < 2 (top left), and in 
the rapidity intervals \y\ < 1 and 1 < \y\ < 2 for the Y(1S) (top right), Y(2S) (botttom left) 
and Y(3S) (bottom right). The uncertainties on the points represent the sum of the statistical 
and systematic uncertainties added in quadrature, excluding the uncertainty on the integrated 
luminosity (11%). 
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Figure 7: (Left) Y(1S) rapidity-differential cross section in the transverse momentum range 
Pj < 30 GeV/c (data points) and normalized PYTHIA prediction (line). The uncertainties on the 
points represent the sum of the statistical and systematic uncertainties added in quadrature, 
excluding the uncertainty on the integrated luminosity (11%). (Right) Y(nS) cross-section ratios 
as a function of pj in the rapidity range \y\ < 2. 
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Figure 8: Differential cross sections of the Y(nS) as a function of pj in the rapidity range \y\ < 2, 
and comparison to the pythia predictions normalized to the measured pj-integrated cross 
sections; Y(1S) (left), Y(2S) (middle), and Y(3S) (right). The pythia curve is used to calculate 
the abscissa of the data points [28 J. The uncertainties on the points represent the sum of the 
statistical and systematic uncertainties added in quadrature, excluding the uncertainty on the 
integrated luminosity (11%). 
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The Y(nS) integrated cross sections are expected to increase with v^. We compare our mea- 
surement of the Y (IS) integrated cross section in the central rapidity region \y\ < 1 to previous 



measurements from the D0 and CDF experiments [4. 5J in Table 11 Previous measurements 
are restricted to the range px < 20 GeV/c and \y\ < 0.4 for CDF and \y\ < 1.8 for D0. Under the 
assumption that the cross section is uniform in rapidity for the measurement range of each ex- 
periment, the cross section we measure at ^/s = 7 TeV is about three times larger than the cross 
section measured at the Tevatron. Although our measurement extends to higher pj than the 
Tevatron measurements, the fraction of the cross section satisfying pj > 20 GeV/ c is less than 
1% and so can be neglected for this comparison. We compare the normalized px-differential 
cross sections at the Tevatron to our measurements in Fig. |9] 
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Figure 9: Comparison of the CMS differential Y(nS) cross sections as a function of pj, normal- 
ized by (JxoT = E(dcr/dpT)Apx/ to previous measurements; Y(1S) (left), Y(2S) (middle), and 
Y(3S) (right). 



9 Summary 

The study of the Y(nS) resonances provides important information on the process of hadropro- 
duction of heavy quarks. In this paper we have presented the first measurement of the Y(nS) 
differential production cross section for proton-proton collisions at a/s = 7 TeV. Integrated 
over the range px < 30GeV/c and \y\ < 2, we find the product of the Y(1S) production 
cross section and dimuon branching fraction to be cr(pp — > Y(IS)X) ■ S(Y(1S) H^h ) = 
7.37 ± 0.131q 42 i 0.81 nb, where the first uncertainty is statistical, the second is systematic, and 
the third is associated with the estimation of the integrated luminosity of the data sample. Un- 
der the assumption that the cross section is uniform in rapidity for the measurement range of 
each experiment, the cross section we measure at = 7 TeV is about three times larger than 
the cross section measured at the Tevatron. The Y(2S) and Y(3S) integrated cross sections and 
the Y(1S), Y(2S), and Y(3S) differential cross sections in transverse momentum in two regions 
of rapidity have also been determined. The differential cross-section measurements have been 
compared to previous measurements and PYTHIA. Finally, the cross section ratios of the three 
Y(nS) have been measured. 

The dominant sources of systematic uncertainty on the cross-section measurement arise from 
the tag-and-probe determination of the efficiencies and from the integrated luminosity nor- 
malization. Both will be reduced with additional data. The cross sections obtained in this work 
assume unpolarized Y(nS) production. Assuming fully-transverse or fully-longitudinal po- 
larization changes the cross section by about 20%. With a larger accumulated data sample, it 
will become possible to perform a simultaneous measurement of the polarization and the cross 
section. This work provides new experimental results which will serve as input to ongoing 
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theoretical investigations of the correct description of bottomonium production. 
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Table 3: The imcorrected Y signal yield, fit quality (normalized 'x^, obtained by comparing the 
fit PDF and the binned data; the number of degrees of freedom is 112), and average weight 
in px intervals for \y\ < 2. The mean of the px distribution in each interval is also given. 
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14-17 


15.4 


l.Z 


708-1-1 8 


U.DU 




17-20 


18.3 


U.O 


lUDztl J 


U.D4 




20-30 


23.3 


U.O 


1 no-i-i '3 

iU7±iO 


U.DU 




sum 






/oZd±1oj 






combined fit 




78n7-l-1 T3 




Y(2bj 


0-2 


1.3 


1 7 
1./ 


^/^8-l-/l 1 
Joozt41 


n /1 7 
U.4/ 




2-4 


2.9 


1.3 


591 ±50 


0.40 




4-6 


4.9 


0.9 


416±40 


0.32 




6-9 


7.3 


1.1 


424±38 


0.33 




9-12 


10.3 


1.1 


257±25 


0.41 




12-16 


13.6 


1.3 


121±16 


0.46 




16-20 


17.7 


1.0 


63±11 


0.55 




20-30 


22.5 


0.8 


39±9 


0.60 




sum 






2279±91 






combined fit 




2270±91 




Y(3S) 


0-3 


1.8 


1.5 


397±51 


0.47 




3-6 


4.3 


1.0 


326±47 


0.37 




6-9 


7.3 


1.1 


264±36 


0.35 




9-14 


11.0 


1.2 


207±25 


0.43 




14-20 


16.3 


1.2 


83±14 


0.52 




20-30 


23.4 


0.8 


49±10 


0.61 




sum 






1324±84 






combined fit 




1318±84 
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9 Summary 



Table 4: The product of the Y(nS) production cross sections, a, and the dimuon branching frac- 
tion, B, measured in pj bins for \y\ < 2, with the assumption of unpolarized production. The 
statistical uncertainty (stat.), the sum of the systematic uncertainties in quadrature (Zsyst.)/ and 
the total uncertainty (A(7; including stat., Ylsyst> ^^id luminosity terms) are quoted as relative 
uncertainties in percent. Values in parentheses denote the negative part of the asymmetric un- 
certainty. The fractional change in percent of the cross section is shown for four polarization 
scenarios: fully-longitudinal (L) and fully-transverse (T) in the helicity (HX) and CoUins-Soper 
(CS) frames. 



Vt 


a-B 


stat. 




^a 


HX-T 


HX-L 


CS-T 


CS-L 


(GeV/c) 


(nb) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 






Y(1S) 




|y| <2 










0-30 


7.37 


1.8 


8(6) 


14(13) 


+16 


-22 


+13 


-16 


0-1 


0.30 


8 


10(7) 


17(15) 


+16 


-22 


+17 


-23 


1-2 


0.90 


5 


9(6) 


15 (14) 


+16 


-20 


+19 


-24 


2-3 


1.04 


5 


8(6) 


14(13) 


+15 


-20 


+19 


-24 


3^ 


0.88 


6 


9(7) 


15 (14) 


+18 


-23 


+18 


-23 


4-5 


0.90 


6 


8(6) 


15 (14) 


+18 


-23 


+16 


-21 


5-6 


0.82 


6 


8(6) 


15 (14) 


+17 


-23 


+13 


-19 


6-7 


0.64 


7 


8(5) 


15 (14) 


+17 


-22 


+11 


-16 


7-8 


0.51 


7 


8(6) 


15 (14) 


+16 


-22 


+7 


-10 


8-9 


0.33 


8 


8(6) 


16 (14) 


+16 


-22 


+4 


-5 


9-10 


0.25 


8 


9(6) 


16 (15) 


+15 


-21 


+2 


-1 


10-12 


0.36 


6 


8(5) 


15(14) 


+15 


-21 


-1 


+3 


12-14 


0.18 


8 


9(5) 


16(14) 


+15 


-20 


-3 


+7 


14-17 


0.14 


9 


10(6) 


17(15) 


+14 


-19 


-4 


+9 


17-20 


0.06 


12 


10(6) 


19(17) 


+13 


-18 


-4 


+10 


20-30 


0.06 


12 


10(6) 


19 (17) 


+12 


-17 


-4 


+10 






Y(2S) 




|y|<2 










0-30 


1.90 


4.2 


9(6) 


15 (13) 


+14 


-19 


+12 


-15 


0-2 


0.25 


12 


11 (9) 


20 (19) 


+14 


-19 


+17 


-22 


1-A 


0.48 


8 


12 (10) 


18(17) 


+12 


-17 


+18 


-23 


4-6 


0.41 


10 


10(8) 


18(17) 


+16 


-22 


+15 


-20 


6-9 


0.41 


9 


10(7) 


17(16) 


+15 


-21 


+9 


-13 


9-12 


0.21 


10 


9(6) 


17(16) 


+14 


-20 


+1 


-0 


12-16 


0.09 


13 


10(7) 


20 (19) 


+14 


-19 


-2 


+6 


16-20 


0.04 


18 


11(8) 


24 (23) 


+12 


-18 


-4 


+9 


20-30 


0.02 


23 


20 (18) 


32 (32) 


+12 


-17 


-5 


+11 






Y(3S) 




\y\ < 2 










0-30 


1.02 


6.7 


11(8) 


17(15) 


+14 


-19 


+10 


-13 


0-3 


0.26 


14 


10(8) 


21 (19) 


+13 


-18 


+16 


-22 


3-6 


0.29 


14 


18 (17) 


26 (25) 


+13 


-18 


+16 


-21 


6-9 


0.24 


14 


11 (8) 


21 (19) 


+15 


-20 


+10 


-13 


9-14 


0.16 


12 


10(8) 


19(18) 


+15 


-20 


-1 


+2 


14-20 


0.05 


17 


11(8) 


23 (22) 


+13 


-18 


-4 


+9 


20-30 


0.03 


20 


12 (9) 


26 (25) 


+11 


-16 


-4 


+9 
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Table 5: The product of the Y(nS) production cross sections, a, and the dimuon branching 
fraction, B, measured in px bins for |y| < 1 and 1 < |y| < 2, with the assumption of impolar- 
ized production. The statistical uncertainty (stat.), the sum of the systematic uncertainties in 
quadrature (Sgyst.), and the total uncertainty (Act; including stat., I^gyst.' luminosity terms) 
are quoted as relative uncertainties in percent. Values in parentheses denote the negative part 
of the asymmetric uncertainty. The fractional change in percent of the cross section is shown for 
four polarization scenarios: fully-longitudinal (L) and fully-transverse (T) in the helicity (HX) 
and CoUins-Soper (CS) frames. 



Pt 


a-B 


stat. 




Act 


HX-T 


HX-L 


CS-T 


CS-L 


(GeV/c) 


(lib) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 




















0-30 


4.03 


1.3 


8 (6) 


14 (12) 


+16 


-22 


+13 


-16 


0-2 


0.70 


5 


9(7) 


15 (14) 


+14 


-19 


+18 


-24 


2-5 


1.54 


4 


10 (9) 


15(15) 


+14 


-20 


+18 


-23 


5-8 


1.02 


5 


7(6) 


14 (13) 


+18 


-23 


+8 


-12 


8-11 


0.44 


6 


7(5) 


15 (14) 


+18 


-23 


-1 


+2 


11-15 


0.23 


7 


8(5) 


15(14) 


+18 


-23 


-4 


+10 


15-30 


0.11 


9 


8(6) 


16 (15) 


+15 


-20 


-5 


+12 




















0-30 


1.03 


2.9 


9 (6) 


15(13) 


+14 


-19 


+12 


-15 


0-3 


0.29 


10 


17(16) 


22 (21) 


+10 


-14 


+17 


-22 


3-7 


0.41 


10 


16(15) 


21 (21) 


+13 


-18 


+14 


-19 


7-11 


0.22 


11 


9(7) 


18 (17) 


+17 


-22 


+1 


-2 


11-15 


0.06 


16 


9(6) 


21 (20) 


+17 


-22 


-4 


+8 


15-30 


0.04 


17 


9(7) 


22(21) 


+14 


-20 


-5 


+11 










\y\ < 1 










0-30 


0.59 


4.8 


11 f8l 


16(15) 


+14 


-19 


+10 


-13 


0-7 


0.38 


11 


25 (24) 


30 (29) 


+11 


-16 


+14 


-19 


7-12 


0.15 


15 


10 (8) 


21 (20) 


+16 


-22 


+1 


-1 


12-30 


0.07 


14 


10 (8) 


20 (20) 


+15 


-21 


-4 


+10 










l< y\<2 








0-30 


3.55 


1.2 


8 (6) 

o yoj 


14(12) 


+16 


-22 


+13 


-16 


0-2 


0.55 


7 


11 (9) 


17(16) 


+18 


-24 


+18 


-23 


2-5 


1.39 


4 


9(7) 


15 (14) 


+20 


-25 


+18 


-23 


5-8 


0.97 


5 


9(5) 


15(13) 


+16 


-22 


+14 


-18 


8-11 


0.37 


7 


10(6) 


16(14) 


+13 


-19 


+6 


-8 


11-15 


0.18 


8 


10 (6) 


17(15) 


+11 


-17 





+1 


15-30 


0.10 


9 


11(6) 


18(16) 


+10 


-16 


-3 


+6 






Y(2S) 




l<\y\<2 








0-30 


0.93 


3.0 


9(6) 


15 (13) 


+14 


-19 


+12 


-15 


0-3 


0.21 


15 


24 (23) 


30 (29) 


+17 


-23 


+17 


-23 


3-7 


0.44 


9 


12 (8) 


18 (17) 


+17 


-22 


+17 


-22 


7-11 


0.19 


12 


11(8) 


20 (18) 


+13 


-18 


+9 


-12 


11-15 


0.06 


17 


11(7) 


23 (21) 


+11 


-17 


+1 





15-30 


0.03 


21 


13(9) 


27 (26) 


+10 


-16 


-3 


+7 






Y(3S) 




1< [i/l <2 








0-30 


0.40 


4.9 


11 (8) 


16(15) 


+14 


-19 


+10 


-13 


0-7 


0.24 


18 


29 (27) 


36 (35) 


+16 


-22 


+17 


-22 


7-12 


0.10 


22 


13 (10) 


28 (27) 


+13 


-18 


+10 


-13 


12-30 


0.06 


17 


11(8) 


23 (22) 


+10 


-15 


-2 


+5 
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9 Summary 



Table 6: The product of the Y(1S) production cross section, a, and the dimuon branching frac- 
tion, B, measured in rapidity bins and integrated over the pj range p ^ < 30 GeV/ c, with the as- 
sumption of unpolarized production. The statistical uncertainty (stat.), the sum of the system- 
atic uncertainties in quadrature (Egyst.), and the total uncertainty (A(7; including stat., Egyst.' ^^id 
luminosity terms) are quoted as relative uncertainties in percent. Values in parentheses denote 
the negative part of the asymmetric uncertainty. The fractional change in percent of the cross 
section is shown for four polarization scenarios: fully-longitudinal (L) and fully-transverse (T) 
in the helicity (HX) and CoUins-Soper (CS) frames. 



Iyl 


a-B 


stat. 


Esyst. 


Act 


HX-T 


HX-L 


CS-T 


CS-L 




(nb) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 






Y(1S) 




Pt < 30GeV/c 








0.0-2.0 


7.61 


1.8 


8(6) 


14 (13) 


+16 


-22 


+13 


-16 


0.0-0.4 


1.62 


3 


8(6) 


14(13) 


+15 


-19 


+13 


-17 


0.4-0.8 


1.52 


4 


9(8) 


15 (14) 


+17 


-22 


+11 


-15 


0.8-1.2 


1.77 


4 


9(7) 


14(13) 


+16 


-22 


+9 


-12 


1.2-1.6 


1.47 


4 


9(7) 


15 (13) 


+17 


-23 


+12 


-16 


1.6-2.0 


1.23 


4 


11(7) 


16 (14) 


+18 


-23 


+20 


-24 
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Table 7: Relative values of the systematic uncertainties on the Y(nS) production cross sections 
times the dimuon branching fraction, in pj intervals for |y| < 2, assuming unpolarized pro- 
duction, in percent. The following abbreviations are used: A, etrig,id/ Sp, Ap^, A^tx, ApsR/ T&P, 
and ej/i/,,Y/ for the systematic uncertainties arising from imperfect knowledge of the acceptance, 
trigger and muon identification efficiencies, momentum scale, the production pj spectrum, the 
efficiency of the vertex quality criterion, the modeling of FSR, the T&P method, and the bias 
from using the ]/xp to determine single-muon efficiencies rather than the Y. The uncertainties 
associated with the background PDF are in the column labeled BG, while the signal PDF, the 
fitter, tracking efficiency, and effects arising from the efficiency binning are combined in the 
colimm labeled add. Values in parentheses denote the negative part of the asymmetric uncer- 
tainty. The luminosity uncertainty of 11% is not included in the table. 



Pt (GeV/c) 


A 


^trig,id 


Sp 


Ap^ 


^vtx 


^FSR 


T&P 




BG 


add. 




Y(1S) 


l3/l<2 


uncertainties are in percent 










0-30 


0.5 (0.5) 


7.5 (4.6) 


0.3 (0.3) 


0.6 


0.7 


0.7 


0.0 


0.9 


0.5 


3.0 


0-1 


0.4 (0.4) 


8.3 (5.4) 


0.1 (0.1) 


0.2 


1.1 


0.8 


0.5 


0.8 


3.4 


3.1 


1-2 


0.4 (0.4) 


7.8 (5.2) 


0.2 (0.2) 


0.6 


0.6 


0.7 


0.2 


1.1 


1.8 


3.0 


2-3 


0.5 (0.5) 


7.3 (4.7) 


0.6 (0.6) 


0.3 


0.3 


0.8 


0.1 


1.1 


1.5 


3.0 


3-4: 


0.6 (0.6) 


7.3 (4.8) 


0.6 (0.6) 


0.1 


0.4 


0.8 


0.0 


1.1 


3.7 


3.0 


4-5 


0.6 (0.6) 


7.4 (4.5) 


0.4 (0.3) 


0.3 


0.7 


0.7 


0.0 


0.9 


2.3 


3.0 


5-6 


0.6 (0.6) 


7.4 (4.3) 


0.2 (0.3) 


0.5 


1.0 


0.7 


0.0 


0.7 


0.5 


3.0 


6-7 


0.6 (0.6) 


7.4 (4.1) 


0.2 (0.3) 


0.7 


1.1 


0.6 


0.1 


0.7 


0.4 


3.0 


7-8 


0.6 (0.6) 


7.7 (4.7) 


0.1 (0.1) 


1.0 


0.7 


0.6 


0.2 


0.8 


1.0 


3.1 


8-9 


0.6 (0.6) 


7.4 (4.2) 


0.0 (0.1) 


1.2 


0.7 


0.5 


0.0 


0.7 


1.0 


3.0 


9-10 


0.5 (0.5) 


7.8 (4.3) 


0.1 (0.0) 


1.3 


0.9 


0.5 


0.2 


0.6 


1.9 


3.1 


10-12 


0.5 (0.5) 


7.4 (3.7) 


0.1 (0.1) 


1.4 


0.8 


0.5 


0.2 


0.6 


0.2 


3.0 


12-14 


0.5 (0.4) 


7.9 (4.0) 


0.2 (0.1) 


1.6 


0.9 


0.5 


0.1 


0.6 


0.3 


3.1 


14-17 


0.4 (0.4) 


8.5 (4.2) 


0.1 (0.1) 


1.6 


0.9 


0.5 


0.3 


0.6 


2.2 


3.1 


17-20 


0.4 (0.4) 


8.9 (4.4) 


0.1 (0.1) 


1.8 


0.8 


0.4 


0.5 


0.7 


0.1 


3.6 


20-30 


0.3 (0.3) 


8.9 (4.3) 


0.1 (0.1) 


1.6 


0.7 


0.5 


0.3 


0.6 


0.1 


3.5 




Y(2S) 


\y <2 


















0-30 


0.6 (0.6) 


8.3 (4.9) 


0.3 (0.3) 


0.7 


0.8 


0.8 


0.0 


1.0 


1.9 


3.2 


0-2 


0.5 (0.5) 


8.3 (5.2) 


0.2 (0.2) 


0.5 


0.6 


0.8 


0.4 


0.6 


6.8 


3.3 


2-A 


0.7(0.7) 


8.3 (5.4) 


0.7(0.8) 


0.2 


0.3 


1.0 


0.1 


1.5 


8.0 


3.3 


4-6 


0.8 (0.7) 


7.9 (4.7) 


0.4 (0.4) 


0.4 


1.1 


0.8 


0.0 


0.9 


5.2 


3.3 


6-9 


0.7(0.7) 


8.6 (4.8) 


0.1 (0.1) 


1.0 


1.2 


0.7 


0.2 


0.9 


1.7 


3.5 


9-12 


0.5 (0.5) 


8.4 (4.2) 


0.1 (0.1) 


1.5 


1.0 


0.5 


0.2 


0.8 


0.9 


3.6 


12-16 


0.4 (0.5) 


8.8 (4.6) 


0.1 (0.1) 


1.6 


0.9 


0.5 


0.3 


0.8 


2.0 


4.0 


16-20 


0.3 (0.4) 


8.3 (4.1) 


0.2 (0.1) 


1.7 


1.0 


0.5 


0.4 


0.5 


0.0 


6.5 


20-30 


0.3 (0.3) 


9.1 (4.4) 


0.1 (0.1) 


1.7 


0.8 


0.5 


0.2 


0.3 


0.0 


17.3 




Y(3S) 


y| <2 


















0-30 


0.7(0.6) 


8.6 (4.7) 


0.3 (0.3) 


0.8 


0.8 


0.8 


0.1 


1.0 


3.4 


5.4 


0-3 


0.5 (0.5) 


8.5 (4.4) 


0.4 (0.5) 


0.5 


0.4 


0.9 


0.2 


0.6 


1.7 


5.7 


3-6 


0.9 (0.8) 


9.1 (5.4) 


0.7(0.7) 


0.3 


0.9 


1.0 


0.0 


1.7 


14.1 


7.3 


6-9 


0.7(0.7) 


8.9 (4.8) 


0.2 (0.2) 


1.1 


1.0 


0.7 


0.0 


1.0 


2.2 


5.6 


9-14 


0.5 (0.5) 


7.5 (4.1) 


0.1 (0.1) 


1.5 


0.8 


0.5 


0.3 


0.7 


0.4 


6.1 


14-20 


0.4 (0.4) 


8.8 (4.5) 


0.2 (0.1) 


1.7 


0.8 


0.5 


0.3 


0.6 


3.4 


5.9 


20-30 


0.3 (0.3) 


8.8 (4.1) 


0.1 (0.1) 


1.6 


0.8 


0.5 


0.5 


0.5 


0.3 


8.3 



24 



9 Summary 



Table 8: Relative values of the systematic uncertainties on the Y(nS) production cross sections 
times the dimuon branching fraction, in pj intervals for |y| < 1 and 1 < |y| < 2, assuming 
unpolarized production, in percent. The following abbreviations are used: A, etrig,id/ Sp, Ap-^, 
^vtx/ ^FSR/ T&P, and £j/^,y/ for the systematic uncertainties arising from imperfect knowledge 
of the acceptance, trigger and muon identification efficiencies, momentum scale, the produc- 
tion Pj spectrum, the efficiency of the vertex quality criterion, the modeling of FSR, the T&P 
method, and the bias from using the ]/ip to determine single-muon efficiencies rather than the 
Y. The uncertainties associated with the background PDF are in the column labeled BG, while 
the signal PDF, the fitter, tracking efficiency, and effects arising from the efficiency binning are 
combined in the column labeled add. Values in parentheses denote the negative part of the 
asymmetric uncertainty. The luminosity uncertainty of 11% is not included in the table. 



Pt (GeV/c) 


A 


^trig,id 


Sp 


Ap^ 


^vtx 


^FSR 


T&P 




BG 


add. 



Y(1S) lyl < 1 uncertainties are in percent 



0-30 


0.5 (0.5) 


7.5 (4.6) 


0.3 


(0.3) 


0.6 


0.7 


0.7 


0.0 


0.9 


0.5 


3.0 


0-2 


0.4 (0.4) 


7.7(5.6) 


0.3 


(0.3) 


0.6 


0.5 


0.7 


0.7 


1.3 


1.6 


3.0 


2-5 


0.6 (0.6) 


7.1 (5.2) 


0.7 


(0.7) 


0.2 


0.1 


0.9 


0.3 


1.5 


6.2 


3.0 


5-8 


0.7(0.7) 


6.5 (4.4) 


0.3 


(0.3) 


0.8 


0.8 


0.7 


0.4 


1.0 


0.3 


3.0 


8-11 


0.5 (0.5) 


6.4 (3.9) 


0.0 


(0.0) 


1.3 


0.5 


0.5 


0.3 


0.7 


0.6 


3.0 


11-15 


0.5 (0.4) 


6.6 (3.8) 


0.1 


(0.1) 


1.5 


0.7 


0.5 


0.3 


0.6 


0.4 


3.0 


15-30 


0.3 (0.4) 


7.1 (4.2) 


0.1 


(0.2) 


1.6 


0.6 


0.5 


0.3 


0.5 


0.9 


3.0 


Y(2S) 


y <i 


0-30 


0.6 (0.6) 


8.3 (4.9) 


0.3 


(0.3) 


0.7 


0.8 


0.8 


0.0 


1.0 


1.9 


3.2 


0-3 


0.6 (0.5) 


8.2 (6.0) 


0.6 


(0.6) 


0.5 


0.1 


1.0 


0.7 


1.4 


13.9 


3.4 


3-7 


0.8 (0.8) 


7.7(5.2) 


0.6 


(0.7) 


0.4 


0.6 


1.0 


0.4 


1.5 


13.1 


3.4 


7-11 


0.6 (0.6) 


7.7 (4.9) 


0.1 


(0.0) 


1.3 


0.8 


0.6 


0.3 


1.0 


1.7 


3.4 


11-15 


0.5 (0.5) 


7.3 (4.4) 


0.1 


(0.1) 


1.6 


0.8 


0.5 


0.3 


0.7 


1.6 


3.6 


15-30 


0.3 (0.4) 


7.4 (4.3) 


0.1 


(0.2) 


1.7 


0.5 


0.5 


0.2 


0.4 


1.9 


4.2 




Y(3S) 


|y|<i 




















0-30 


0.7(0.6) 


8.6 (4.7) 


0.3 


(0.3) 


0.8 


0.8 


0.8 


0.1 


1.0 


3.4 


5.4 


0-7 


0.8 (0.8) 


8.8 (5.9) 


0.6 


(0.7) 


0.5 


0.5 


1.0 


0.5 


1.7 


22.9 


5.4 


7-12 


0.6 (0.6) 


7.6 (5.0) 


0.0 


(0.0) 


1.4 


0.6 


0.6 


0.1 


0.9 


2.2 


5.7 


12-30 


0.4 (0.4) 


7.1 (4.0) 


0.2 


(0.1) 


1.6 


0.7 


0.5 


0.3 


0.5 


0.2 


6.0 




Y(1S) 


l<|y|<2 


















0-30 


0.5 (0.5) 


7.5 (4.6) 


0.3 


(0.3) 


0.6 


0.7 


0.7 


0.0 


0.9 


0.5 


3.0 


0-2 


0.4 (0.4) 


8.2 (4.6) 


0.0 


(0.1) 


0.3 


1.1 


0.6 


0.2 


0.6 


7.3 


3.0 


2-5 


0.5 (0.5) 


7.7 (4.0) 


0.3 


(0.3) 


0.2 


0.9 


0.7 


0.3 


0.5 


4.3 


3.0 


5-8 


0.6 (0.6) 


8.4 (4.2) 


0.1 


(0.2) 


0.7 


1.2 


0.6 


0.5 


0.6 


0.4 


3.0 


8-11 


0.6 (0.5) 


8.9 (4.4) 


0.0 


(0.1) 


1.3 


1.1 


0.5 


0.6 


0.7 


1.7 


3.1 


11-15 


0.4 (0.5) 


9.1 (4.2) 


0.1 


(0.2) 


1.6 


0.9 


0.5 


0.8 


0.6 


0.1 


3.0 


15-30 


0.4 (0.5) 


10.6 (4.3) 


0.2 


(0.3) 


1.7 


1.1 


0.5 


0.9 


0.8 


2.0 


3.1 



Y(2S) 



1 < 



1/1 < 2 



0-30 


0.6 (0.6) 


8.3 


(4.9) 


0.3 


(0.3) 


0.7 


0.8 


0.8 


0.0 


1.0 


1.9 


3.3 


0-3 


0.5 (0.5) 


7.8 


(3.8) 


0.2 


(0.2) 


0.3 


1.2 


0.7 


0.2 


0.5 


21.9 


3.8 


3-7 


0.7(0.7) 


9.4 


(4.9) 


0.3 


(0.3) 


0.4 


1.3 


0.8 


0.4 


0.6 


5.4 


3.4 


7-11 


0.6 (0.6) 


9.6 


(4.8) 


0.0 


(0.0) 


1.2 


1.1 


0.6 


0.7 


0.8 


4.4 


3.5 


11-15 


0.5 (0.5) 


9.7 


(4.8) 


0.1 


(0.2) 


1.6 


0.9 


0.5 


1.1 


0.7 


1.5 


4.6 


15-30 


0.4 (0.4) 


9.5 


(3.8) 


0.2 


(0.2) 


1.7 


1.3 


0.4 


1.0 


0.7 


3.4 


7.2 




Y(3S) 


1 < 


|y|<2 


















0-30 


0.7(0.6) 


8.6 


(4.7) 


0.3 


(0.3) 


0.8 


0.8 


0.8 


0.1 


1.0 


3.4 


5.5 


0-7 


0.7(0.6) 


8.6 


(2.9) 


0.4 


(0.4) 


0.4 


1.1 


0.8 


0.5 


0.3 


26.5 


6.5 


7-12 


0.7(0.6) 


9.3 


(4.1) 


0.0 


(0.0) 


1.3 


1.0 


0.6 


1.1 


0.6 


6.2 


6.7 


12-30 


0.4 (0.4) 


9.4 


(4.3) 


0.1 


(0.2) 


1.7 


1.1 


0.5 


0.9 


0.5 


1.5 


5.9 
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Table 9: Relative values of the systematic uncertainties on the Y(1S) production cross section 
times the dimuon branching fraction, in rapidity intervals for pj < 30GeV/c, assuming unpo- 
larized production, in percent. The following abbreviations are used: A, etrig,id/ Sp, Ap^, Avtx/ 
^FSR/ T&P, and ej/^,Y/ for the systematic uncertainties arising from imperfect knowledge of the 
acceptance, trigger and muon identification efficiencies, momentum scale, the production pj 
spectrum, the efficiency of the vertex quality criterion, the modeling of FSR, the T&P method, 
and the bias from using the ]/tp to determine single-muon efficiencies rather than the Y. The un- 
certainties associated with the background PDF are in the column labeled BG, while the signal 
PDF, the fitter, tracking efficiency, and effects arising from the efficiency binning are combined 
in the column labeled add. Values in parentheses denote the negative part of the asymmetric 
uncertainty. The luminosity uncertainty of 11% is not included in the table. 



\y\ 


A 


^trig,id 


Sp 


Ap^ 


^vtx 


^FSR 


t&p 




BG 


add. 




Y(1S) pT<30GeV/c 


uncertainties are in percent 






0.0-2.0 


0.5 


7.5 (4.6) 


0.3 (0.3) 


0.6 


0.7 


0.7 


0.0 


0.9 


0.5 


3.0 


0.0-0.4 


0.6 


6.8 (4.9) 


0.4 (0.4) 


0.7 


0.3 


0.7 


0.6 


1.5 


0.1 


3.0 


0.4-0.8 


0.6 


6.8 (4.7) 


0.4 (0.4) 


0.6 


0.3 


0.7 


0.3 


1.1 


5.4 


3.0 


0.8-1.2 


0.5 


7.5 (4.9) 


0.3 (0.3) 


0.6 


1.0 


0.7 


0.1 


0.7 


2.9 


3.0 


1.2-1.6 


0.5 


7.7 (4.0) 


0.2 (0.2) 


0.6 


1.2 


0.6 


0.2 


0.5 


4.0 


3.0 


1.6-2.0 


0.6 


9.3 (4.0) 


0.0 (0.1) 


0.6 


0.9 


0.6 


0.9 


0.6 


5.0 


3.0 



Table 10: The ratios of Y(nS) cross sections for different Y pj ranges in the unpolarized scenario. 
The first uncertainty is statistical and the second is systematic. The ratios are independent of 
the luminosity normalization and its uncertainty. 



Pt (GeV/c) Y(3S) / Y(1S) Y(2S) / Y(1S) 



0-30 


0.14 ±0.01 ±0.02 


0.26 ± 0.02 ± 0.04 


0-3 


0.11 ±0.02 ±0.02 


0.22 ± 0.03 ± 0.04 


3-6 


0.11 ±0.02 ±0.03 


0.25 ± 0.03 ± 0.05 


6-9 


0.17 ±0.03 ±0.03 


0.28 ± 0.04 ± 0.04 


9-14 


0.20 ± 0.03 ± 0.03 


0.33 ± 0.04 ± 0.05 


14-20 


0.26 ± 0.07 ± 0.04 


0.35 ± 0.08 ± 0.05 


20-30 


0.44 ±0.16 ±0.08 


0.36 ± 0.14 ± 0.06 



Table 11: Y(1S) cross-section measurements at several center-of-mass collision energies. The 
first uncertainty is statistical, the second is systematic, and the third is associated with the 
luminosity determination. 



Exp. 


V's 


|(pp ^Y(IS)X) 


rapidity 




(TeV) 




range 


CDF 


1.8 


0.680 ± 0.015 ± 0.018 ± 0.026 nb El 


\y\ < 0.4 


D0 


1.96 


0.628 ± 0.016 ± 0.065 ± 0.038 nb |H 


\y\ < 0.6 


CMS 


7.0 


2.02 ± 0.03;g {^ ± 0.22 nb (this work) 


\y\ < 1-0 
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